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It was recently argued that the observed PVLAS anomaly can be explained by chameleon held 
I theories in which large deviations from Newton's law can be avoided. Here we present the predictions 

>— ^ , for the dichroism and the birefringence induced in the vacuum by a magnetic held in these models. 

We show that chameleon particles behave very differently from standard axion-like particles (ALPs). 
{^JQ' We find that, unlike ALPs, the chameleon particles are confined within the experimental set-up. As 

^ , a consequence, the birefringence is always bigger than the dichroism in PVLAS-type experiments. 
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PACS numbers: 14.80.-j, 12.20.Fv 



I. INTRODUCTION 

Light scalar fields are common in theories of physics beyond the standard model. These scalar fields can couple 
to both the standard model fields as well as to new types of matter. Since experiments have not yet detected new 
1^ I forces, it means that the force mediated by these scalar fields is either very weak (with a strength less than gravity) 
. or short-ranged in the laboratory. In string theory, for example, there are many moduli fields which couple to matter 
with gravitational strength. The chameleon mechanism provides a way to suppress the forces mediated by these scalar 
. fields via non-linear field self- interactions and interactions with the ambient matter [J, As a result, the masses 
^ ' of the scalar fields become dependent on the ambient matter density; this is the reason why these fields have been 
^-H dubbed chameleon fields. If the observed accelerated expansion of the universe is due to a chameleon-like field, it has 

■ interesting cosmological consequences [1, 0| . 

I An obvious way to look for chameleon fields are gravitational experiments in different environments P, . Addi- 
. . tionally, it was recently pointed out that chameleon fields are a natural way to reconcile the PVLAS and CAST results 
t ' 0. In 2006, it was reported that the PVLAS experiment had detected light polarization rotation in the vacuum in 
, the presence of a magnetic field, j6j] . Recently there has been a lot of activity concerning the theoretical interpretation 

■ of this detection, see e.g. [3|-[23| for recent work. In particular, this finding could be seen as evidence for the presence 
of an axion-like particle (ALP) or milli-charged particles. The PVLAS 2006 results can be explained if the mass of the 
ALP is of order mALP ~ 10~^ eV and the inverse coupling constant to two photons is M « 10^ GeV. These results 
are in direct conflict with the CAST results, because, for these parameters, axions emitted in the sun should have 
been detected by this experiment. However, as was pointed out in [5|, assuming the values for m and M given above, 
both results can be explained if the particle behaves like a chameleon field, because the mass of the field inside the 
PVLAS experiment would then be very different from the mass of the field inside the Sun. In the Surr, the chameleon 
mass is so high that chameleon particles cannot be photo-produced. Additionally, it was shown in I5| that the new 
force mediated by the chameleon field is not in conffict with current experiments. 

Although in a context different from chameleon theories, the potential resolution of conflict between the CAST and 
the PVLAS experiments provided by particles with an environmentally dependent mass was also discussed in Ref. 
Q. It was also pointed out that if the particle mass outside the cavity was much larger than inside, then ALPs would 
be reflected in the same way as the photons and not escape from the Fabry-Perot cavity. As an implication, there 
would be no signal found in "light shining through the wall" experiments iS]. As the BMV collaboration seems to 
have observed [9] . It should be noted that in Ref. Q the ALPs were assumed to reflect in exactly the same way as the 
photons; we see below that this is not the case if the ALP is a chameleon. The predictions of the chameleon model 
are therefore very different from the ones of reflecting ALP model consider in Ref. ^ . 

Recently new PVLAS results have been reported These new results do not confirm the previously reported 
rotation. Indeed they find no evidence for either rotation or ellipticity with a 2.3 T magnetic field. The initial PVLAS 
results were found using a 5.5 T magnetic field. With a 5.5 T field, the new results additionally show no evidence 
for any rotation, however a non-zero ellipticity is still detected. In the context of standard ALPs (but not necessarily 
chameleon fields), however, the detected ellipticity is excluded with a 99% confidence by the null result found with 
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the 2.3 T magnetic field. Although it might be an artefact, the ellipticity for B = 5.5 T can be explained with a 
chameleon model with m « 10^"^ eV and M « 10^ GeV. Future experiments will give us decisive clues later this year. 

In this paper we study the behaviour of chameleon particles inside optical experiments similar to PVLAS. The 
aim of this paper is to derive the expression for the rotation and the ellipticity of the laser light. As we will show, 
chameleon particles and standard ALPs behave very differently. Firstly, because the mass of the particles depend 
on the environment and they do not have enough energy, chameleon particles cannot leave the experimental region. 
In contrast to ALPs, which do leave the interaction region, chameleon particles are therefore trapped inside the 
experimental set-up. Secondly, because the mass of the particles varies inside the experiment, chameleon particles 
and photons reflect differently. As we will discuss, this results in a very different form for the expressions for the 
predicted rotation and ellipticity. 

The paper is organized as follows: In Section |TT] we present the chameleon model in more detail. In Section IIIII 
we discuss how the chameleon mass behaves inside the experimental set-up. The behaviour of chameleon particles in 
experiments like PVLAS is described in Section [TVl The predictions made by the chameleon model are discussed in 
Section |Vl We conclude in Section IvH Mathematical details can be found in the Appendices. 



II. THE CHAMELEON MODEL 



Chameleon theories are essentially scalar field theories with a self-interaction potential and couplings to matter; 
they are specified by the action 



1 p'P/M 



^4 

+ Y^S^ie^'^'^g,,,^^) (1) 

where the Sm and ij^m are respectively the matter actions and matter fields. The couplings to matter and the 
electromagnetic sector are specified by the mass scales Mi and M respectively. For simplicity we take the coupling 
of the scalar field to matter and to the electromagnetism to be universal (i.e. Mi = M). In more general theories 
one would not expect a truly universal coupling. Our conclusions are, however, not affected by this assumption. The 
coupling to matter implies that particle masses in the Einstein frame depend on the value of 

m[c^) = e-^/^^mo (2) 

where mo is the bare mass as appearing in Sm- The strength of the chameleon to matter coupling is given by 



where Mpi — I/vSttG w 2.4 x 10^^ GeV. As we show in Section ITVl if a theory such as this is to be detected by axion 
searches one must require that M ^ Mpi. This implies that on scales smaller than Sc/m^, where is the mass of 
(/), the chameleon force between matter particles is 2{Mpi/M)^ ^ 1 times stronger than their mutual gravitational 
attraction. If the mass, m^, of is a constant then one would then have to require that ^ ImeV otherwise 
such a theory would already be ruled out by experimental tests of Newton's law. Another, and potentially far more 
interesting prospect, is that the mass of the scalar field grows with the background density of matter. In high density 
regions it could then be large enough to satisfy the constraints coming from tests of gravity, whilst being small enough 
to produce detectable alterations to the standard physical laws in low density regions. Scalar fields that have this 
property are said to be Chameleon fields. In addition to the couplings to matter, chameleon fields have non-linear 
self- interactions described by a potential V{(j)). Assuming an exponential coupling to matter of the form given by Eq. 
P]), a scalar field theory will have a chameleon mechanism, for some range of (j), provided that: 

V'{<j)) < 0, V" > 0, V"'{(f>) < 0, (4) 

where V — dV/dcj). Whether or not the chameleon mechanism is strong enough to evade current experimental 
constraints depends partially on the details of the theory, i.e. 1^(0) and M, and partially on the initial conditions 0]. 
For exponential couplings and a potential of the form 

A 4+n 

y (</>) = A^ eM^Vr) « A^ + — - (5) 
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the chameleon mechanism can in principle hide the field such that there is no conflict with current laboratory, solar 
system or cosmological experiments [l|, y] ■ Importantly, the chameleon mechanism is strong enough in such theories 
to allow strongly coupled theories with M <C Mpi to have remained thus far undetected 

The first term in the potential, V , corresponds to an effective cosmological constant whilst the second term is a 
Ratra-Peebles inverse power law potential. If one assumes that (f) is also responsible for late-time acceleration of the 
universe then one must require A « 2.3 x 10^^^ GeV. 

The evolution of the chameleon field in the presence of both matter and an external magnetic field is determined 
by the effective potential: 

V,s{<jy)^V{^)+p,se^'^ (6) 

where 

Poff = Pmattcr + (7) 

and B is the magnetic field. As a result, even though V has a runaway form, the effective potential has a minimum 

at = (?!)min(Poff) where 

Kff (0mi„) = = V^'(0,„i„) + — f Pmattcr + ^ j • (8) 

In the presence of ambient matter and a magnetic field, the field evolves towards this minimum. The mass of small 
perturbations in about ^min is given by 

n+4\ 1/2 

(9) 



/ A"+4\ 

TO0((^min) = (n(n+l)-— 



and we have that 

2nA4+"M 



2/3mattor + B 



1/(1+") 

(10) 



In some circumstances, however, (j) is unable to change quickly enough to actually reach 4>miu- In particular, this 
behaviour can occur inside a low density cavity. If the radius, R, of the cavity is too small, then does not reach its 
effective minimum and instead, as we show in Section Hill '-^ 2/R for R <^ 2/m^{(l)min)- Since laboratory searches 
for vacuum magnetic dichroism and birefringence generally employ such a cavity, one must be particularly wary of 
this behaviour when making predictions for what such experiments should detect. In particular the dependence of 
on B and Pmattcr inside the experiment depends on the details of the set-up. 
The chameleon mass, m^, depends therefore on a number of factors. In searches for dichroism and birefringence in 
a vacuum, the mass depends on the details of the experimental set-up itself: the size of the cavity, the magnetic field 
B and the density of matter in the laboratory vacuum. The mass is not a fundamental parameter; the fundamental 
parameters of our model are A, M and n. It is also important to note that, not only is not a fundamental 
parameter, but that it is generally very different in different parts of experiments. We discuss the behaviour of in 
laboratory searches for axion-like-particles (ALPs) in the next section. 



III. THE CHAMELEON MASS IN THE LABORATORY 

In laboratory searches, such as the PVLAS @, and Q&A 24l| experiments, for the dichroism and birefringence 
induced by a magnetic field, light propagates in a Fabry-Perot cavity with radius R. The interaction region, i.e. the 
region where the magnetic field is turned on (_B 7^ 0), has length L. To increase the strength of any signal the light 
is refiected N times, and the mirrors are located a distance d from either end of the interaction region. We label 
the density of the vacuum matter inside the cavity by Pgas- For example in the PVLAS experiment L — 100cm, 
d ~ 270cm, R = 12.5mm and pgas ~ 2 x 10~''gcm~^. Before we can make predictions for the dichroism and 
birefringence we need to know the value of along the path of the photon. Although we are primarily concerned 
with power-law type potentials, in most of the discussion in this section we do not assume any particular form of V, 
we only require that it satisfies the chameleon field theory conditions given by Eq. 

The cavity is a cylinder with radius R. Outside of the cavity (i.e. in the walls of the cavity and the surrounding 
magnet), (j) must lie close to its effective minimum, which we label by = If this were not the case then the walls 
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of two such cavities would feel a force that would be 2{M pi / M)"^ ^ 1 times their gravitational attraction. Such a 
force is easily ruled out by experimental tests of gravity. We define r to be the radial distance from the centre of the 
cavity and define 0o = (t'i^ = 0). We approximate the potential inside and outside the cavity as a quadratic function. 
Outside the cavity, the expansion is around the minimum where VJ,'ff (0oo) — 0. Inside we expand around <f>^ which is 
left unknown. This leads to two equations for the field <f). For r > R we have 

</)" + i0'-m^(0-(/.^) = O, (11) 



whereas for r < R the equation reads 

1 



+ ^(/.'~m2(0-</,o)-Kffo, (12) 



where l/'^g = V"'(0o) + Pcse'^"/^^ ■ Note that in most cases 0o is not (/)c where V^'ff(0c) = 0. We define = V'{(j)c) 
and Vq = V'{(f)o), we can then write V^f^ o ~ ~ K'- 

The non-singular solution for r < ^ is a combination of the Bcsscl functions Jq and Nq. However, since A^o diverges 
logarithmically at the origin, we ignore the term containing A'o and therefore the solution inside the cavity reads 

cf,^CMimor)+(bo- ^° ■ (13) 

TTIq 

On the other hand, for r > R, the solution is 

(j) = A{Jo{imoor) ~ iNo{inioor)) + 0oo (14) 

Matching both solutions and their first derivatives at r = i? gives the following conditions for A and C: 

mnJ^{imoR){(j)oo - (f>o + jVf'^ - K0/™o) 

mocJoiimoR){Jlj - iN^){imooR) - moJ'a{imaoR){Jo - iNo){imaoR) ' 

and 

^ _ nioo (Jq - iN^){im^R) 
Too jQ{imoR) 

Now, since (j) must lie very close the effective minimum inside the cavity walls, we must have 

niooR > 1- 

The solution inside the cavity then simplifies drastically 

, , , v'~v' 



Jo[imoR) 

Evaluating this equation at r = and imposing (/)(r — 0) — (po leads to 



Joiimor) + (jjo - ° 2 ^ (15) 



iMimoR) - 1) (16) 



TO, 







There are now two relevant situations: 



A. rricR > 1 

Since toq < TOc, moR ^ 1 in this case. It follows from Eq. that Vq w and therefore that toq w TOc 
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B. rricR < 1 

If Vq/V^ ^ 1 ^ ^(1) that would require mg ~ rric and by the matching condition moR « rricR ^ C'(l)- If nicR <C 1 
we must therefore have Vq/V^ — 1 ^ 1, and so 

1 + Ja[nnaR) = -— . 

I ''o I 

The right hand side of this equation is generaUy 0(1) or smaller and so mg-R = C(l) which implies mg <C rrioc and 
ipoo 00 due to the runaway form of the potential. For the inverse power-law potential Vq/toq = ~(f>Q/{n + 1), and 
the matching equation therefore gives 

Jo{imoR)^n + 2, (17) 

which implies, as expected, that mgi? = C(l) for n — 0(1). Expanding the Bessel function for small mgi? <C 4 in the 
last equation gives 

moR~2VWTT. (18) 

A slightly better approximation for 0(1) values of n is given by: 

TOoi?~ 2%/2(\/^rr2- 1)^/^ (19) 

Note that if rricR ^ 0(1) then we could not ignore the term in the matching condition, but since < we would 
have: 

J{imoR) - 1 > ?i + 2, 

and so mgi? is always larger than the value defined by Eq. (|17p . In summary: we have found that there are two 
relevant cases for the chameleon mass. In the first case, the length scale set by the chameleon mass inside the cavity 
(1/mo) is much smaller than the size of the experiment. In the second case, the length scale I/toq is 0(i?), and we 
found an approximate relation Eq. (fl7|) between the mass inside the cavity, R and the properties of the potential, 
encoded here in the power n. Eq. ([TT]) also defines the smallest possible value of mg. It is clear that for 0(1) values 
of n + 1 we cannot have mg <^1/ R. 

IV. CHAMELEON OPTICS IN A CAVITY 

It was shown in Ref. that to avoid constraints on solar ALP production from the CAST experiment, one 
must require that in bodies with densities of the order of lOgcm"^ the chameleon mass satisfies > 10** eV. In 
experiments such as BRFT, PVLAS, and Q&A the photon beam typically has a frequency lo = 0(l)eV. In the 
mirrors then 3> uj and so the chameleon field cannot escape the cavity. Indeed the propagation of the chameleon 
field outside the cavity is exponentially attenuated implying that the cavity mirrors also act as perfect mirrors for the 
chameleon field. This is at odds with the usual assumption that ALPs escape from the cavity without any reflection. 
The standard expressions for the rotation (dichroism) and the ellipticity (birefringence) must therefore be modified. 

In this section we study the propagation of a beam of light in the presence of a chameleon field and derive expressions 
for the predicted rotation and ellipticity. For simplicity we initially assume that there is no distance between the end 
of the interaction region and the mirrors, and that the chameleon and photon fields reflect in the same way. Whilst 
neither of these assumptions are generally true, and there are important effects associated with the violation of each 
of them, the calculation is much simpler and easy to follow if we make them. We say more about what occurs when 
these assumptions are dropped in subsection IIV Dl below. 

We assume that the scalar field mixes with the orthogonal polarization to the magnetic field. This system is a 
two- level system with two states \P > and \S > for the photon and the scalar in the absence of the magnetic field. 
For the system of differential equations satisfied by the photon and the scalar, see the Appendices. When a magnetic 
field is turned on, the two states mix and the eigenstates are 



1+ > = cos^lS* > +isin6'|P > 
I- > = cose\P > +isme\S > 
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where 



tan 26* : 



2Buj 
Mm? 



The eigenvalues for the above system are given by 



For small 9 we get 



and 



In particular we find 



and 



which is crucial in the following. 



7 2 2 2 

k<=bj — m 



cos 26* ± 1 
2 cos 26* 



fc+ = W - 



2w 



2/32 



m^9 
2uj 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



A. Free propagation 

We consider first the situation where the electromagnetic wave and the chameleon propagate freely. Assuming that 
the state at a given origin is |P > (z = 0), then state at a further position z is given by 



\P > [z] = cos A:_2: cos 6'|— > — i cos k^z sin 9\+ > 
This mixed state can be expressed in terms of the free scalar and photon as 

\P > (z) = (cos^ 9cosk^z + sir? 9 cos k^z)\P > +isin0cos0(cosfc_z — cosfc_|-z)|S' > 
The photon part for small 9 and using the expansion of k± (|24l25p . is given by 



(26) 



^(z) 



1 - 29^ sin^ 



cos UJZ + 



2lj 



2uj 



from which we identity the attenuation and the phase shift 



a = 29^ sin^ 



4a; 



(5 = 



m 



2/j2 



2u) 



-z — a sm ■ 



2uj 



(27) 



(28) 



^From those, dividing by two, one gets the rotation and the ellipticity for an incoming laser beam with a 45 degree 
polarization. We finish the subsection defining a quantity which will be very useful below 



-^coh 



2uj 



(29) 



measuring the coherence of the system. 
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B. Propagation in a cavity 

The only difference between the free propagation case and a propagation within a cavity of size L is that we need 

to sum over periodic copies shifted by nL . The mirrors are not perfect, so only N < oo passes actually occur. 
Mathematically, this system is equivalent to the coupled wave equations in a periodic box. 
Let us define first 

a„(z) = 26»2 sin^ ^—^ i (30) 

4cj 



and 



S„(z) = — — {z + nL) - sm (31) 



2uj " ' 2uj 

The photon wave function is then given by summing 

N-l 

fl}~[ 

(^ — n.„ \ cos fc ( 7. -A- n.T:\ -\- (fl sin 

2u 



iP[z) = ^ (1 - a„) cos k- {z + nL) + 6*^ sin sin fc_ {z + nL) (32) 



n=0 

In a perfect cavity, the waves have a resonance , i.e. wL = 2'Kp. Therefore 



tp{z) = ^ (1 - a„) cos ( + ^ — {z + nL)\ 

n=0 ^ ^ 



2 . m'^{z + nL) 



+9 sin sin k- {z + nL) 

2lj 



Using the small 9 approximation, one obtains 

N-l 

N 



2uj 2uj 

n=0 

Notice this is nothing but a standing wave 

JV-l \ / N-l 



N ^ ^ ' \ N 

\ n=0 / \ n=0 

So the phase shift and the attenuation at the end of the cavity z = L are given by 

N-l , Af-1 



For larger values and using sin a; < a;, we have the inequality St < Sn- 
Similarly the attenuation is 



(33) 



il){z) = ( 1 - ^ X! ) cosw-z (34) 

^V^2 . m^{z + nL) m^9\ . 
+ 22(^ sm - -^{z + nL) ) smuz 



n=0 n=0 

These sums can be exactly computed, which give 



' {N + 1) L siniV^+sin^-sin(7V + l)^ \ 

2 2coh 2Af(l-cos^) I ^ ' 



„/ l + cosAr-^-cos-^-cos(iV + l)-^\ 
I 2Ar(l-cos-f-) i ^ ' 
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C. Coherence 



The previous expressions can be easily interpreted in terms of the coherence of the photon-scalar system. We focus 
on the case where the coherence length and the length of the experiment are commensurate. Define the number of 
coherent passes as PL = 2TTZcoh- Using this one finds that 

a„+p = a„, d„+p = On ^ — - — PL (39) 

The attenuation is then given by 

N/P-l p-i N/P-1 p-i 

1 ^ 



^ E E«"+j^ = ^ E E«" (40) 

^■=0 n=0 j=0 n=0 



and therefore 

p-i 

1 



1 

P E (41) 



p 

n=0 

Hence the attenuation depends only on the coherence length, after P passes, the waves are not coherent any more. 
Using the previous formulae we find 

ar = 6*2 (42) 

and for the rotation 

, ^ 

rotation/pass = (43) 

which was first derived in [2^. 

In order to get the ellipticity one notes first that the phase shift is given by 

JV/P-lP-l W/P-lP-l 2/12 \ 

- ^ E E - ^ E E W + ^-^^^j (44) 

j=0 n=0 j=0 n=0 ^ ' 

and therefore 

ST = PLi--l)^HP+l)0^^^ (45) 



Using the definition of P one finds 



ST = n[—-l\6' + 2n^-e' (46) 



For large P this simplifies to 

^rr^ — IT I 

^P 

This can also be written for N/P 1 as 

NLrn^ 



5T^^(^-l]e' (47) 



5t = — 0^ (48) 



Which again was found in [26|. The ellipticity per pass is now 



ehipticity/pass ^ — ^ ) (49) 



and one finds that the ellipticity per rotation is 



rotation \P J \ 27rZcoh 

In particular, for large N, the ellipticity is always much larger than the rotation. This fact is still true when the 
non-interacting zone of length d is taken into account. It is a crucial prediction of chameleon theories. 
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D. Phase Shifts 



The above calculation has been performed under the assumptions that there is no gap between the ends of the 
interaction region and the mirrors and that the chameleon field reflects off the mirror in precisely the same way as 
the photon does. Unfortunately, neither of these assumptions are generally true. 

In the PVLAS, BRFT and Q&A experiments, there is always some gap between the end of the interaction region 
and the mirrors. In PVLAS this gap is 270 cm long. Outside the interaction region both the chameleon field and the 
photon field propagate freely, however since the chameleon field is massive it travels more slowly than the photon. As 
we show in Appendix [Bl this leads to the chameleon field picking up a phase shift relative to the photon field 
when the field returns to the interaction region. This phase shift alters the formulae for the rotation and ellipticity. 
We find that if there is a distance d between the mirror outside the interaction region, then Am ~ mid/uj. In addition 



Very fast (step-like) change in ttie Chameieon Mass More realistic - 0(1 )/d ctiange in Chameleon Mass 




> 

distance from surface of mirror distance from surface of mirror 



FIG. 1: Illustration of the difference between a sharp (step-like) change in the chameleon mass at the surface of the mirror and 
the more realistic ~ 0{l/d), for d S> l/rric, behaviour, where d is the distance from the surface of the mirror. In both of 
these sketches, rric is the chameleon mass inside the mirror and rUb is the chameleon mass far from the mirror. The dotted line 
indicates the surface of the mirror. 

to the phase shift Am , there is generally an additional phase shift A,, due to the manner in which the chameleon field 
reflects off the mirror. Consider a mirror placed at x = 0. If behaved liked a step function, see Figure [H where 
7710 = m< for a; < and = m> for x > and m> <ti m<, then the wave function of a chameleon wave with 
wave-number k incident on the mirror would, in x > 0, have the form: 

6(j) ~ ScjjQ sin(fc2:). 

When X < 0, (j) decays away exponentially. However, in realistic theories (j) is continuous across x — and so is 
too. Indeed, in Appendix [Bl we show, via a similar calculation to the one performed in Section IIIII above, that for 
l/"^< ^ ^ ^ l/'7i> we have « 0{l)/x (see Figure [T|). 

It should be noted that this behaviour is not specific to the inverse power-law potentials considered here, but is 
generic to almost all chameleon theories (see Appendix |B] for a complete discussion) . A scalar wave with frequency lu 
switches from oscillatory to exponentially decaying behaviour when = w^. In chameleon theories this transition 
therefore occurs at a distance Xr = 0{1/uj) from the mirror. Chameleon waves therefore behave as if they reflect not 
at a; = but sd x = Xr > 0- Schematically this alters the form of the chameleon wave function in x > Xr to: 

6(1) w 6(1)0 sin(fc(a; — Xr)). 

As is shown in Appendix |Bl this early reflection produces another phase shift. A,., in the chameleon wave relative to 
the photon wave when they return to the interaction region. We find that for inverse square law potentials: 



where n > 0. 
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The presence of a total phase shift A « + A^, alters the formulae for the rotation, Ac^a and ellipticity, ip and 
these are calculated in Appendix [X] Both of the mentioned effects will exist, to some degree in all theories in which 
the ALP does not escape the cavity. We find that when A 7^ and provided that: 



Buj sin 



in (^m'^L/Au 



Mml sin (a/2 + mlL/Auj] 



«1, 



(51) 



we predict 



Aip 
sin 2ip 





(52) 



and 



sm2ip 




(53) 



where 



Ga(x) 



and 



The phase A is given by 



2sin(A/2)sin(x/2) 
sin(A/2 + a;/2) 
sin2(a;/2) 
sin2(A/2 + a;/2)' 

_ sin((iV + l)x) 
~ {N + l)sinx ' 



(54) 
(55) 

(56) 

(57) 
(58) 

When A = 0, Ga = and Ha — 1- We always have |(5Ar(a;)| < 1, and for N large 5n{x) is strongly peaked about 
X = rrnr for any integer m and (5jv(a;) ^ 1 otherwise. In many situations one finds that m?^{L/2 + d)/2uj ^ 1 and 
tan(Ar/2) 0(1). In these circumstances, the expressions for the ellipticity and rotation simply greatly: 



mid 



where A^. depends on the potential V{<j}). For inverse power-law potential (n > 0), A^ is given by 

■nn 



Aip 
siniy9 



32M2sin2(%)' 



sin^ 16M2tan(%^)' 

which are both independent of m^. In this limit, the ratio of the rotation to the ellipticity is: 

Alp _ 1 
~ ~ 7Vsin(Ar) ■ 



(59) 
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V. PREDICTIONS 

Having described the physics and derived the basic formulae for rotation and elhpticity induced by the coupUng of 
a chameleon field to the electromagnetic sector, we now study the predictions of the chameleon model. To simplify 
the analysis, we have assumed that the chameleon couples to all matter types (including photons) with the same 
strength. In this case, there is a lower bound on the energy scale M coming from the contribution of the chameleon 
to the anomalous magnetic moment of the muon and electron. The contribution is of order (mg.^/M)^ WA - order 
for the chameleon contribution to the anomalous magnetic moment of the muon to be small enough, M has to be 
bigger than M k, 10^ GeV. We will therefore concentrate here on the case M > 10* GeV. We should mention here 
that if the couplings differ from species to species, this constraint on the coupling to the electromagnetic sector could 
be relaxed. 

The predictions for rotation and ellipticity are shown in Figure [21 One feature of the chameleon model is that the 
elhpticity is generally predicted to be much larger than the rotation. This can be viewed as a generic prediction of 
chameleon theories. In general, a large ellipticity to rotation ratio occurs (for large N) whenever: 

• The ALPs are reflected by the mirrors rather escaping through them. 

• The ALPs from previous passes are not coherent with the photon field. 

If the ALPs were to escape then each pass would give the same contribution to the rotation and ellipticity, these add 
up and thus both Aip and would be proportional to A^. However, if the first of the above conditions is satisfied 
then this is not the case, since ALPs from previous passes interact with the photon field and thus the way in which 
the rotation and ellipticity build up is altered. If, additionally, ALPs from previous passes are not coherent with 
the photon field, then a large degree of cancellation of both the ellipticity and rotation occurs. For large N, the 
cancellation in the rotation, Aip, is almost exact, and as result Atp is almost independent of N. 

Cancellation of contributions to the ellipticity also occurs. However, there is always at least one contribution which 
does not cancel but builds up as the number of passes increases. This contribution results from the fact that the 
component of the photon field that interacts with the ALP propagates more slowly in the interaction region than 
its non- interacting component. As a result the relative phase of the interacting and non- interacting components is 
shifted. This contribution to the ellipticity is unaffected by the decoherent ALPs. There is, therefore, always at 
least one contribution to the ellipticity that is proportional to N. Thus the combination of reflecting ALPs and their 
becoming decoherent results in the ellipticity but not the rotation growing as TV. The rotation to ellipticity ratio is 
therefore ~ 1/N ^ 1 (see Eq. ([55| for a better estimate). 

If the ALP in question is a chameleon field then both reflection and decoherence almost certainly occur. Reflection 
occurs because the mass of the chameleon inside the mirror, TOc, is much larger than the mass far outside, mb, and the 
energy of beam, w, used in experiments is almost always <C rric. Decoherence occurs because, as we found in Appendix 
IHl at a distance x {l/rric <C a; ^ l/Tib) from the surface of the mirror, the chameleon mass behaves as ^ 0{l)/x. The 
reflection of the chameleon does not, therefore, occur on the surface of the mirror but a distance Xr ~ 0{l)/uj from 
the mirror's surface. This, we found, resulted in an 0(1) phase shift between the reflected chameleon and photon 
fields. Generally then, the chameleon fields from previous passes are both reflected and are not coherent with the 
photon field. Whilst, one could probably construct a chameleon theory where one or both of the above conditions 
did not hold, it would not be particularly generic and might also have problems satisfying gravitational and other 
laboratory constraints on such theories. 

It is therefore a generic feature of chameleon theories that the rotation to ellipticity ratio is 0{1/N) ^ 1. This 
means that, according to the chameleon model, it is easier to detect the ellipticity than the rotation. 

In Section IlIII we showed how, for a given experimental set-up and choice of V{(l)), the chameleon mass, m^, in 
the interaction could be calculated from the chameleon-to-matter coupling if A and n are known. We found that 
there were two regimes. If the density of the vacuum, Pgas, and the radius, R, of the cavity were small enough and 
M large enough then ^ 0{2)/R for small n. Alternatively, if M is small enough or pgas / R large enough, then 
the chameleon mass depends on (pgas + B'^ /2)/M. The chameleon mass and therefore the predicted ellipticity and 
rotation are therefore highly dependent on the set-up of the axion search experiment. In the chameleon model, there 
is no reason to expect that two different experiments should detect the same ellipticity and rotation. 

A. Predictions for PVLAS, Q&A, BMV and BRFT 

As we argued in Section [Hi we expect A « 2.3 x 10^'^ eV if the energy density of the chameleon field is to be 
associated with dark energy. Ideally, one would determine both A and n from an experimental detection. However, 
since the status of the only (i.e. PVLAS) detection reported to date is unclear, we now consider the specific predictions 
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FIG. 2: Predictions for rotation (upper plot) and ellipticity (lower plot) in the chameleon model as a function of M and m^. 
The PVLAS set-up is used {L = 100 cm, d — 270 cm, u = 1.2eV, B — 5 T and ip = tt/4). Furthermore we have chosen n = 1. 

for the chameleon model for A w 2.3 x 10"^ eV and n = 1. The vacuum used in the PVLAS set-up is very good, with 
a density of about 2 x lO^^'' gcm^'^, whereas the density of the Q&A experiment's vacuum is significantly higher: 
8.5 X lO^'^gcm^^; as a result, for M ~ 10'' - 10^° GeV, the chameleon mass in the PVLAS experiment ~ 0{2)/R 
and hence independent of M. In the Q&A experiment, however, is, for given M, both larger than it is in the 
PVLAS set-up and depends on {pgas + B^/2)/M. The larger value of in the Q&A set-up means that the ellipticity 
and rotation predicted by the chameleon model are far smaller than those predicted for the PVLAS experiment. The 
chameleon model predictions for both the PVLAS and Q&A set-up are shown in Figure [31 The thin dotted lines 
in these plots show the 2.3 T PVLAS 2007 [lO| and Q&A 95% confidence limits on the rotation and ellipticity. For 
M > lO'' GeV it is only the 2.3 T PVLAS upper bound on the ellipticity that provides a useful constraint on M: we 
must require M > 2 x 10® GeV if ri = 1 and A«2.3xl0^'^ eV. A similar limit of M is found for other values of n. 

The PVLAS experiment performs better than Q&A as a probe for chameleon fields because of the high quality 
vacuum it employs. The upcoming BMV experiment uses a similar high quality vacuum with pressure < lO^^mbar 
[l^H^I- This experiment will additionally use a higher strength magnetic field than PVLAS, employ a greater number 
of passes and additionally promises a higher precision. Figure 2] shows the predicted rotation and ellipticity signals 
for PVLAS (with both B = 2.3 T and 5.5 T) and BMV again with n = 1 and A = 2.3 x lO^^gy. We see that BMV 
should be able to detect, or rule out, such chameleon theories with Af < 3 x 10^ GeV which represents an order of 
magnitude improvement over PVLAS. 

For completeness, we show in Figure [5] the predicted rotation and ellipticity signals for PVLAS {B = 2.3 T) and 
BRFT (for two different number of passes A^) again with n ~ 1 and A = 2.3 x 10"'^ eV. Notice that, although 
chameleon models predict a higher rotation and ellipticity within the BRFT set up, we see that BRFT should only 
be able to detect, or rule out, such chameleon theories with M < 8 x 10^ GeV. This is worse than that provided by 
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FIG. 3: Predictions for rotation (left) and ellipticity (right) in the chameleon model as a function of M for A = 2.3 x 10"^ eV 
and n = 1. Predictions for the 2.3 T PVLAS (L = 100 cm, d = 270 cm, uj = 1.2eV, B = 2.3T, pg^s = 2 x lO^^gcm"^ and 
ip = 7r/4) and Q&A (L — 60 cm,d — 145 cm, uj = 1.2eV, B = 2.3T, pgas ~ 8.5 x 10~^gcm~^ and (p — tt/4) set-ups are shown. 
The thin-dotted lines show the 95% confidence upper bounds on both the rotation and the ellipticity. 




FIG. 4: Predictions for rotation (left) and ellipticity (right) in the chameleon model as a function of M for A = 2.3 x 10"^ eV 
and n = 1. Predictions for the B = 2.3 T and B = 5.5 T PVLAS (L = 100 cm, d = 270 cm, uj = 1.2 eV, Pgas = 2 x 10""gcm"^ 
and ip — 7r/4) and BMV {L = 50 cm,d = 85 cm, u) = 1.2 eV, B = 11.5 T, pgas ~ 10~^''gcm~^ and ip = ■7i"/4) set-ups are shown. 
The thin-dotted lines show the 95% confidence upper bounds on both the rotation and the ellipticity. 



PVLAS. Even though the BRFT and PVLAS set-ups use vacuums of similar quality, the smaller number of passes in 
the BRFT experiment result in it placing a much weaker bound on the ellipticity per pass than the PVLAS bound. 



B. Dependence of predictions on n 

In Figures [2]|4] we took n — I. Theories with different values of n do, however, lead to different predictions. Figure 
ini shows how the ellipticity and rotation predicted for the PVLAS set-up depend on n. We found above that if 
X = ml{d + L/2)/2lj < 1 and tan(Ar/2) - 0(1) then: 

^ - ^ 16Mnan(l!(.)/2) (^ iV sin(A) A,, 

which depends on n only through Ar{n) — TTn/{n + 2). Thus for values of M such that X ^ 1, which for PVLAS 
roughly corresponds to M ^ 10^ GeV, both the ellipticity and rotation are therefore only very slightly dependent on 
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FIG. 5: Predictions for rotation (left) and ellipticity (right) in the chameleon model as a function of M for A = 2.3 x 10"^ eV 
and n = 1. Predictions for the B = 2.3 T PVLAS [L = 100 cm, d = 270 cm, u = 1.2 eV, pgas = 2 x lO^^^'gcm"^ and tp = tv/4) 
and BRFT (B = 2 T, L = 800 cm, d = 345 cm, lo = 2.41 eV, Pgas ~ 10"^*gcm"^ and (p = 7r/4) set-ups are shown. The 
thin-dotted lines show the 95% confidence upper bounds on both the rotation and the ellipticity. 




Chameleon to matter coupling: M (GeV) Chameleon to matter coupling: M (GeV) 



FIG. 6: Predictions for rotation (left) and ellipticity (right) in the chameleon model as a function of M for A = 2.3 x 10~^ eV 
and different values of n. The PVLAS set-up is used (L = 100 cm, d = 270 cm, lo = 1.2eV, B = 5 T, pgas = 2 x 10~"gcm"^ 
and ip = 7r/4). 



n and scale as /AP. Different 0(1) values of n do not, therefore, alter the magnitude or M dependence of either 
tp or for M > 10'' GeV. If M < 0(10'^) GeV however, then X > 0(1) then both ■(/' and Aip depend strongly on 
which in turn depends strongly on Pgas, B and n for such values of M. We can clearly see this transition from 
strong to weak n dependence in Figure [SI 



C. Has a chameleon field already been detected? 



The PVLAS 2007 results found no evidence for any ellipticity with B = 2.3 T. However, at B = 5.5 T a non-zero 
ellipticity was detected: ip = (9.7 ± 1.3) x 10^* with 45000 passes |l(|. This is equivalent to an ellipticity per pass 
of: (2.2 ± 0.3) X 10~^^. If, as is the case for standard ALPs, ip (x B'^ then such an ellipticity at 5.5 T implies that at 
2.3 T, PVLAS should have found ip = (1.7 ± 0.2)10"^ which is, in fact, ruled out with 99% confidence. A B"^ scahng 
of ip therefore implies that the signal detected with B = 5.5 T must be of instrumental origin (lo| . 

In chameleon theories, however, nr^ can depend on B and so a B^ scaling of the ellipticity is not assured. We 
noted above that ii X = m?,{2d + L)/2uj <C 1 then both tp and scale as B^ . If we wish to reconcile the 5.5 T 
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detection with the 2.3 T null result we must therefore concentrate on C(l) or greater values of X , which for PVLAS 
corresponds to > 3 x 10~^eV ^ 2/i? « 3 x 10~^ eV. We are therefore well into the region where depends 
strongly on both pgas + /2. Now grows with _B^, but it is generally the case that decreases as to^ grows, 
and therefore ^j: generally grows more slowly than for X ~ 0(1). However, we need a faster than B^ growth if 
we are to find a theory that predicts the signal seen aX B — 5.5 T without violating the B — 2.3 T upper bound. 
Fortunately, the necessarily faster than B^ growth in \'4}\ does occur in chameleon field theories when: 

A, ml{B)2d + L 
- + 

for some integer m. Both tp and Ac/) are strongly peaked about values of that satisfy the above equation. In these 
cases all the effects that usually result in the chameleon and photon fields becoming decoherent cancel each other 
out, which results in a greatly amplified signal. For these values of the experiment could be thought of as being 
resonant. The positions of these resonances are highly dependent on the set-up of the experiment. If m^[B — 5.5 T) 
for S = 5.5 T just so happens lies close to such a resonance but ■m^{B = 2.3 T) does not, then we would have 
\^{B — 5.5 T)| ^ (5.5/2.3)^|?/'(i3 = 2.3 T)|. Additionally, since i/' changes sign as one passes through such a resonant 
point, we can ensure that 'ijj{B ~ 5.5 T) > 0. In the context of standard ALPs, V > implies the existence of 
a pseudo-scalar, in the chameleon model, however, this is not the case and ■0 can be both positive and negative 
depending on the value of to^. 

We find that for A = 2.3 x 10"^ eV, there exist values of M for which V'(5.5 T) = 9.7 ± 1.3 x 10"** and |i/'(2.3 T)| < 
1.4 X 10~* for n > 2. In all cases M ~ 0(10^ GeV) and ~ O(10~^ eV) with larger values of n corresponding to 
both larger values of M and smaller values of |V'(2.3 T)|. For example: If n = 2 then M = (1.13 ±0.03) x 10® GeV gives 
V'(5.5T) = 9.7± 1.3 X 10"**, V(2.3T) = -(1.32±0.05) x 10"^, Av5(5.5T) = -(1.4±0.2) x lO-^^ and A(p{2.3T) = 
-(1.6 ± 0.1) X 10"^^. If this is the case then the BMV experiment would measure %p{BMV) w -7.3 x 10"^. If 
however n = 3 then M = (1.58 ± 0.04) x 10® GeV is required and |'(/'(2.3T)| < 4 x 10"^; BMV would measure 
ip{BMV) « —3.5 X 10"^. We have checked that these models would are not ruled out by the BRFT experiment. 
In all cases, the ellipticity and rotation predicted by these models in the context of that set-up is well below the 
BRFT upper bounds. This is due in part to the long length of the BRFT interaction region (8 m) which results in 
its sensitivity being peaked for particles with mass smaller than 1 meV. In contrast, all of the chameleon models that 
reproduce the signal observed by PVLAS have m > 1 meV i.e. outside of the region where BRFT works best.. 

In the chameleon model then -0(5. 5T) = 9.7 ± 1.3 x 10~* is not necessarily excluded by |i/'(2.3T)| < 1.4 x 10"^, 
however it does require that chameleon mass for B = 5.5 T lies very close to a resonant point. Since the position of 
these resonances is highly dependent on the set-up of the experiment this would be a remarkable coincidence. This 
issue will certainly be settled by the upcoming BMV experiment. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have studied how chameleon fields that couple to the electromagnetic sector alter the propagation 
of photons in a vacuum. Just as other axion like particles (ALPs) do, such chameleon fields would induce both 
dichroism (rotation) and birefringence (ellipticity) in a photon beam travelling through a magnetic field. Both of 
these effects could be detected by laboratory searches for ALPs such as PVLAS [e| , Q& A [2!] and BMV f29|, [s^ . The 
mass of a chameleon field depends on its environment; specifically it is larger in backgrounds where the ambient matter 
density is high than it is in those where the background density is low. The mass and coupling strength of standard 
(i.e. non-chameleonic) ALPs are strongly constrained by limits on solar axion production. As was pointed out in 
however, the density dependence of the chameleon's mass implies that in the Sun is generally much larger than 
the value of in the laboratory vacuum. Solar axion production therefore represents a far less stringent constraint 
on chameleon ALPs than it does on standard ALPs. This opens to the door to the prospect that chameleons fields, 
if they exist, may well be detected first by ongoing and upcoming laboratory axion searches. 

The original motivation for this work was the reported detection of light polarization rotation in the vacuum in the 
presence of a magnetic field by the PVLAS experiment Q. If one wishes to explain this detection by the presence 
of a standard ALPs then it must have mass ttialp ~ 1 meV and photon coupling M w 10® GeV. A standard ALP 
with these properties is, however, strongly ruled out by bounds on solar axion production 5], and if the ALP is a 
scalar field then it is additionally ruled out by short-range laboratory tests of gravity. However, as was shown in 0], 
a chameleon field with these properties is not ruled out. If chameleon fields behaved the same as standard ALPs they 
could, at first sight, explain the PVLAS detection. 

In this paper we have shown that, within the confines of experiments like PVLAS, chameleon fields and standard 
ALPs behave very differently. As a photon beam moves through the experiment it generates ALPs. In the absence 
of a chameleon mechanism these ALPs pass through the mirrors at either end of the Fabry-Perot cavity and escape 
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the experiment. Chameleon fields, however, do not escape. As a direct result of their chameleonic properties, the 
chameleon mass in the mirrors is many orders of magnitude larger than in the cavity. The mirrors therefore 
act as a potential well for the chameleon particles and, as was shown above, these particles do not have energy to 
propagate through the mirrors. Not only then do the photons reflect off the mirrors, but so also do the chameleon 
particles. Furthermore, whereas the photons reflect off the surface of the mirror, a beam of chameleon particles with 
energy oj was found to reflect at a distance 0{\/ljj) away from the surface of the mirror. A coherent beam of photons 
and chameleons incident on the mirror is therefore decoherent after reflection. The presence of these two effects means 
that the standard expressions for the rotation and ellipticity no longer hold. In Section [IVI we therefore derived new 
expressions for the dichroism and birefringence induced in a photon beam by the presence of a chameleon field. The 
new expressions for the rotation and ellipticity are given respectively by Eqs. ((5^ and ([55)1 . 

The combination of both the refiection and decoherence of the chameleon particles was seen to lead to the magnitude 
of the predicted ellipticity always being much larger than that of the predicted rotation. Specifically, if the photon 
beam makes N passes through the Fabry-Perot cavity, then the ratio of the rotation to the ellipticity is 0{1/N) <C 1. 
The large ellipticity to rotation ratio is a generic feature of chameleon theories, and so experimental searches for 
birefringence are better placed to detect / rule out the existence of chameleon particles than measurements of rotation. 

In Section |V] we used our expressions for the rotation and ellipticity to make predictions for the P VLAS [6| , Q& A 
and upcoming BMV [2l,[30| experiments. The magnitude of the potentially detectable signal was found to depend 
heavily on the density of the laboratory vacuum. The relatively high density vacuum used in the Q&A experiment 
was seen to result in a much smaller predicted rotation and ellipticity than that for the PVLAS and BMV set-ups. 
If we assume that the non-zero ellipticity detected by PVLAS with B — 5.5 T is instrumental then the most recent 
PVLAS results constrain M > 2 x 10'' GeV for n = 1 and A = 2.3 x 10"^ eV, with a similar bound for other 0(1) 
values of n. BMV will be able to detect or rule out the presence of chameleon fields with A = 2.3 x 10~^eV and 
M ^ 3 X 10^ GeV. On the other hand, if the B = 5.5 T ellipticity signal is physical, then chameleon theories can 
provide an explanation with a mass w 10"^ eV and an inverse coupling M « 10^ GeV. In this case, a large 
ellipticity should be observed by the BMV experiment, providing a decisive tests of chameleon theories. 

In conclusion, chameleon theories can be tested with light propagating in vacuum through a magnetic field as 
long as the coupling to the electromagnetic sector is large enough {M ^ Mpi). Chameleon theories are then an 
additional interesting class of models to be probed by experiments of PVLAS-type. We found that the predictions 
of chameleon theories are substantially different from those of standard ALPs and as such could potentially reconcile 
astrophysical with local tests. These experiments provide additional constraints on the parameter of the theory, which 
are complementary to gravitational and/or Casimir-experiments. 



APPENDIX A: GENERALIZED CALCULATION OF PROPAGATION OF LIGHT AND THE 

CHAMELEON IN A CAVITY 

In Section IIVI we calculated how light and the chameleon field propagate in a cavity under the assumptions there 
is no gap between the mirrors and the ends of the interaction region and if the photon and chameleon field reflect 
in the same manner. In general neither of these two assumptions hold. The calculation performed in Section IIVI is 
very direct and elegant however we found that it becomes significantly more complicated when the assumptions are 
dropped. In the limit where both of the aforementioned assumptions hold, the results found using this alternative 
approach are entirely equivalent to those found in Section [IVI In this appendix our method is based on that used in 

Waves in the photon field. A, and in the chameleon field, (j>, obey: 

□A = Xiii^, (Al) 

n^-ml^ - (A2) 

We take B = Bbx, A — a\\ej; + a±ey. The waves then travel in the z-direction. We take the interaction region, where 
B 7^ 0, to have length L. There is a distance d between the ends of the interaction region and the mirrors. 
For right moving waves cx exp(ifcz), in the interaction region, we have: 

d^a\\ — — fc^O||, 

and 



- dlv = Ur-v, (A3) 
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where 

V — 

Ur = 

and where x = The eigenvectors of Ur are: 



X 



fc2 -kB/M \ 
-kB/M P J ' 



cos 9 \ I — sin 9 

sin 6* / ' \ cos 6* 



with eigenvalues and respectively where: 

■ = -I- m"'- 

2 COS 26* 



2 ,2 2COS2^± 1 _ ^, 



As in Section IrVl we have defined: 



tan 26* 



Mm2 ' 



Eq. (jASP can then be written as: 



where 



cos u — sm t 
sin 9 cos 6* 

and fl — diag(a;^, w^). A similar equation holds for the left-moving modes: 
where 



Q- 



~ cos t/ sm t 
sin 9 cos ( 



In the interaction region then the evolution of the right moving modes is given by: 

v(i)=e— V+(i)v(0), 

and the evolution of the left moving modes by: 

v(t)-e— V_(t)v(0), 

where V± = Q^diag (e-*('^--'^)*, e-*('^+-")*) g±. For t = L we have: 



where: 



c± = cos2 6le-*("=f-'^)^ +sin2 6'e-'('^±-'^)^, (A6) 
&= )i_e-'(-+-M (A7) 



Thus a wave that enters the magnetic field region at z = as vq will, upon exiting, the field region at z = L have 
evolved into: 

v(L) =e-*"^t/+(L)vo. 
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The wave then exits the magnetic field, travels a distance d to the mirror, is reflected, travels a further distance d 
and re-enters the magnetic field at z — L. Outside the interaction region the chameleon field travels more slowly than 
the photon, and furthermore it may reflect at a different point. In general then the chameleon will return to z = L 
having been shifted by a phase relative to the photon field. We define this phase to be A. The assumption made in 
Section ITVl are equivalent to taking A = 0. We calculate this A in Appendix [Bl 

When the fields re-enters the field region, this time moving leftwards, they are therefore (up to an overall irrelevant 
phase factor) in a state: 

v = e-^'^^i?T/+(L)vo, 

where R = diag(l, — exp(— iA)). The — sign is in the last component of R due to the fact that chameleons are scalar 
fields with positive parity, whereas the photon fields are vectors with negative parity. The fields then travels through 
the interaction region, this time moving leftwards, exit the region and reflect back again, so that by the time they 
have returned to z = moving rightwards they are in a state: 

v.oW = e-2-i {RV^RV+)^^o■ 

If A is taken to have an imaginary part then this analysis also allows us to account for any scalar fields that escape 
the cavity. For instance if, as is often taken to be case when the ALP is not a chameleon, the scalar field docs not 
reflect then A = —zoo, and R = diag(l, 0). We deflne: 

F = RV.RV+ = ( ^Ja^ g Ja^ ] , (AS) 



where 



Then after N passes through the interaction region, the fields will be in a state: 

vjv — e I' vq. 



We define 



and find that: 



Zn 



An+i = uAn + iZn, (A9) 
Zn+i = e-'^7v4Ar + e-'^pZN , (AlO) 



and so: 



An+i = {a + Pe~'^) An + e~'^ [7^ - a0\ A^-i- (All) 
We take initial conditions Aq = sin if and Zq = which implies that Ai — a sin ip. We then find the solution: 

An _ f a - a - fi+ ^ 



sm (fi — /i- fi+ — II- 

where 



_ ^ (A12) 



fi+ = a + e-'^^Sii, (A13) 
1^- = e-^^ iP ~ 6ii) , (A14) 



where 
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and 



If A = then: 



We take ^ <C 1 and so: 
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2b 



{a - /3e-»^) c+ - e-'^c- 



li± = e 



We define 5±/2w = a;± - w and (5 = 5+ - For small 9, 5_ = 0{e^ 



iA 



(X < 1. 



After N round trips, a photon which was initial in a state {cosip, sini^)^ has evolved into: 



Oil W \ - { ^(^) + ^"pW) 
a±{N) J ~ {r]{N)sm{ip + Aip{N))e'p'^^'> 



cos 

tan(y+Ay(Af)) ^io(Af) ^.^ ^ 
tan (fi 



where 



tan(v? + A(^(iV))ev(^) An 



tany> 



sin^' 



(A15) 



Since the photon is trapped in the cavity, what one actually measures is a superposition of states which have each 
completed a different number of passes. One therefore measures: 



Measured — , -i 2-^^\'^> — I tan(y+Ay„) -q 



where 



tan(y + Aip^)e'P"' 
tant^ 

We calculate A(j)m and pm to 0{6'^). To this order 

w e-2'^e-'''+^/'^ f 1 + 26'^ 

We define 



1 ^ A 

N + 1 ^ sine' 

k=0 



(l-e-'^)(l-e-">'-^/2")' 

]^ _ Q-iA-i6L/2u; 

■(l-e-*^)(l-e-*''^/2'^) 



]^ _ g—iA—iSL/2u 



(A16) 
(A17) 



ifA(a;) 



sin(a;/2) 



sin(A/2 + a;/2) 



Assuming that 9'^Ha{5L/2u}) <^ 1, we perform the sum over the Ai~ and find that we find: 



tan((/J + A(prn)t 

tanip 



(A18) 



where we have defined 



2sin(A/2)sin(x/2) 
GaW = r 



sin(A/2 + x/2) ' 
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and 



sm{{N + l)x 
{N + l)sina;' 



If A = then Ga = and Ha = 1- The rotation is given by Aipm and the elhpticity by -0 = — pmSin2(^/2. For 
general real A we therefore find: 



and the ellipticity is 



sin 2ip 




2 / A^A Nm^^L 
"2~ ^ 4w 



sin2(y5 



1 



Buj 



2 I Mm! 



2uj 



f NmlL\ 
sin iVA + ^U. 



iVGA 
2uj , 



2w 



<5Ar A + 



2uj 



As a consistency check we can also consider the case where the scalar fields escape after every pass. This is given 
by taking exp(— lA) 0, and in this limit Ha and GAix) sin(a;)2j sin^(a;/2) and so we recover the standard 
formulae: 



escape 



sin 2(p 



Blu 
Mml 



sm 



Aid 



and 



sin 2lp 



1 / Buj 

2 \ Mml 



NmlL ( mlL' 

- iVsin 



2cj 



2uj 



In Appendix [B] we consider the propagation and reflection of the chameleon field outside the interaction region and 
evaluate A. 

APPENDIX B: REFLECTION OF THE CHAMELEON FIELD 

We consider how chameleon waves with frequency lo reflect off a flat mirror that is placed at z — zq, and propagate 
relative to a photon field outside the interaction region. We write (j) — (j)(){z)+d(f>{z, t), where Scj) is a small perturbation 
with frequency w, and (j)o is the background value of (j). We begin by showing that ~ l/{z — zq) near the mirror. 

1. Behaviour of the chameleon mass near a mirror 



We define p.^ to be the density of the mirror and define (f>rn, by: 

We define mm = 'Ti0(0,„). We take the mirror to lie in the region z < zq with its surface at z = zq. Since M <^ Mpi, 
4>Q must be « (pm deep inside the mirror for consistency with experimental tests of gravity. Inside the mirror, (p^ 
obeys 
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and outside the mirror (z > zq), where the density of matter is Pgas7 we have 
We define 0c by: 



PgeLS 



M ' 

and rric = m^{(j)c). Integrating the above equations and assuming that deep inside the mirror (j) — > and d(j)/Az 
we find: 

\(^)^ = V{Mz))-V[c^^)+^-^{Uz)-<P^) z<zo, (Bl) 

\{^)^ ^ V{M^))^V[c^,) + P^{M^)-c^,) z>zo. (B2) 
Matching these equations sd z = zq and using p,„ 2> Pgas we find that at z = zq: 

Mz = w 0„ . (B3) 

^From Eq. (jB2p . we see that outside the mirror, in the region where (0c)/V^'(0o)| 1, we have: 

2 



We now normahze the potential V so that as pgas — > 0, V^((/'c(Pgas)) — 0, i.e. we neglect any constant term in V. 
\V'{(j)c)/V'{(j)o)\ < 1 then implies than V{(j)o) > V{(t>c) and so: 

djl/mp) 

— — ^ 1' (^4) 

where 

^/2(F"(0o))^/' 



(-F"'(</>o))(t/(0o))i/2 



> 0. 



If F oc " then — ^2n(n + l)/(ri + 2)^. We define a0(zo) = and use the shorthand a0((/)(z)) = a0(z), we then 
have: 

1 1 z-zo 1 r*^^-' 1 , 

— = — + — -^ + — -/ — da^. B5 

mo ms 00 (z) a0(z) J^^ toq 

In many theories, the potential is such that changes only very slowly i.e. \d\i\a^l dx\ <ti |c?lnmo/da;| for lo < toq ^ 
m™; when this is true we have: 

«0(^) .T^«^ 

mo ~ ■ — , (B6) 

z - Zo + a^[z) ms 



where a^{z) is slowly varying compared to mo and rris = yJV" {(I)q{z — zq)) is the mass of the chameleon on the surface 
of the mirror. 

Since we have assumed that \V' {(pc) /V {4't))\ ^ 1 then we must have mo(z) 3> mc which implies z — zq ^ i/rnc- 
For l/rrim ^ l/m^ <C Az = z — zo ^ ^/ttIc, Eq. (jB6[) gives mo oc 1/Az. 

This behaviour will occur in any chameleon theory where varies slowly with (j) compared to m^. More generally, 
we have: 

- a0(z) 
Too > 



z - Zo + a^{z)/ms ' 
if 

2V""V" V"V' 

and mo < a0/(z — zo + a^/rris) if fo^ > 1- The assumption that is slowly varying compared to 1/too is equivalent 
to 3|1 — 6^1 ^ 1 for mo < m^ < ms < m^. 
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2. Reflection of chameleon waves 

We now consider the reflections of chameleon waves: Scf) = (50(z)e~*'^*. These evolve according to: 




For most sensible choices of potential we have TOs ^ 0{mm — and we assume that mm S> w; which must 

certainly be the case if w ^ 0{1) eV and the constraints on solar axion production are satisfied [H]- For z < zq then 
we have 64> ~ Cexp(r7ie//(z)(z — zq)) where meff{z) ~ 0{m^{(j3m)) and C is a constant. 

The behaviour of mo for mo S> nic = iTi^{(j>c) is given by Eq. (jB6[) . In the experiments that we consider mc ^ cj. 
The behaviour of mo for mo <C to only effects this reflection calculation at sub-leading order and so, to a first 
approximation we have taken mo(2;) to be given by Eq. (jB6p . Defining x = z — zq + a^/rris, we find that near the 
mirror we have: 

dx'^ \ x"^ 
which has solutions: 

— \fZjX {ciJa{LOx) + C2Na{Ldx)) , 

where Ja and Na are Bessel functions and a — iy'l + 4a^. Near z = zo we have ujx « a^uj/nis = S. Generally 
6 — a^uj/ms ~ 0{uj/mm) <C 1. Matching at z = zq we find: 

C = VSiciM5) + c2N^{S)), (B7) 
mi<5C = ^VS{c,{J^{S) + 2SJ'^{6))+C2{N^{6) + 26K{6))), (B8) 

where mi '--^ 0(ms) is given by mi = dexp(—meff{z)(z — ZQ))/dz\^^^^. The precise value of mi is generally not 
important. We find then that: 



where 

m 




In general S <ti 1 which implies that B{S) <C 1 and so C2 ^ ci. A distance z 3> l/m^{(j)c) from the mirror then S(j) is 
therefore given by: 



-Cisin(^(z-zo)- J)) • (B9) 



We write: 



A chameleon wave that leaves the interaction region, reflects off the mirror is therefore, on its return to the interaction 
region, shifted by a phase A = 2(fc — + relative to a photon that has made the same journey. From Eq. (jB9p 
we have 



A, = I (2a -1) = 1(^1 + 4^-1). (BIO) 

If Ar ^ 1 then the additional phase-shift due to the the chameleon travelling more slowly than the speed of light 
(the 2(fc — Lo)d term in A) might be important. Taking k = y^cij^~~m|(0^, we find: 

mid 



2{uj - k)d ' 
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We therefore we find: 

A = A,, + « I (y^l + 4a2 - 1 



UJ 



Throughout this work we have been primarily concerned with inverse power law potentials: 1^ — oc 0"" for n > 0. 
Valid chameleon theories exist with n — —2m, where to is a positive integer. In a more general scenario the potential 
might, around the point where « lj look locally like cq + ci(/>^", for some co and ci, or any n. For such theories 
if n < — 2 or n > —2/3 then we have: 



A, = 



whereas if —2 < n < —2/3 wc find: 



7r(n + 1) 
"2(?i + 2)' 
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